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We present a theory of phonon- assisted optical two- 
magnon absorption in two-leg spin-ladders. Based on the 
strong intra-rung-coupling limit we show that collective ex- 
citations of total spin S = 0, 1 and 2 exist outside of the 
two-magnon continuum. It is demonstrated that the singlet 
collective state has a clear signature in the optical spectrum. 

71.27.-ha, 75.10.Jm, 75.50.Ee, 75.30.Et, 78.30.Hv 
ITP preprint: NSF-ITP-99121 

Recent design of materials with localized spin- 1/2 mo- 
ments arranged in ladder, as well as frustrated and 
dimerized chain geometries has enormously intensified 
research on low-dimensional quantum magnetism pi. 
Many of the novel materials are likely candidates for 
a spin-liquid ground state and display a gapped spin- 
excitation spectrum. On the two-leg ladder the spin-gap 
is a phenomenon well established by various theoreti- 
cal approaches, i.e. exact diagonalization PJq], strong 
coupling expansion [yj5|, density-matrix renormalization 
group [|j, and bosonization |^. Optical spectroscopy, 
and in particular Raman scattering, has proven to be a 
valuable tool to investigate the spin excitations of the 
new quantum magnets [p[. Apart from Raman scat- 
tering, optical absorption has has been studied, how- 
ever, only on the sm^fe-chain cuprate Sr2Cu03 pjl0| 
where it has provided direct evidence for the gapless two- 
spinon spectrum of the spin-1/2 Heisenberg chain pi[ . 
In Sr2Cu03 absorption is due to phonon assisted two- 
magnon emission (PME) |13,n3l. While several of the 
new spin-ladder materials, eg. SrCu203 and LaCu02.5, 
may allow for PME as well, corresponding studies are 
lacking. In this work we propose a simple theory of PME 
for spin- ladders. We show that bound states in the spin- 
gap which have no analogy in the Heisenberg chain have 
a substantial effect on the optical absorption spectrum. 

We start our discussion with the two-leg spin-ladder 
Hamiltonian 



H = y\s^^iS^i 
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where S°i with a — x,y, z is a spin-1/2 operator on site 
I of leg /x and H is measured in units of J± with A = 
J\\/ J±- In this paper we focus on the limit of large intra- 
rung coupling A — > in which a reformulation of ffl) in 
terms of the rung-spin S" — S^i + 5^; and the operator 



T,° 



'11 



S21 is appropriate 
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/ , S"S^^i - Hi + H2 + H3 
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where Si — 0(1) for a rung singlet (triplet). The rung-spin 
eigenbasis is given by \s) = (| ti> - I iT>)/\/2, \tx) = 
-(I TT) - I ii))/%/2, \ty) = z(| tt) + I ii»/V2, and 
\tz) = (I Ti) + I iT))/'s/2 where the first(second) entry in 
the kets refers to a site on leg T(2)' of the ladder. As 
usual S"\s) — and -S'^ji^) = i£af3j\tf) with the Levi- 
Civita symbol Safs-y- The action of T" on the rung basis 
is given by T"|s) = |i„) and r"|t^) = S<,f^\s). 

For vanishing inter-rung coupling the ground state 
of (||) is a pure rung-singlet product-state |) = <<S)i\si). 
The excited states are products of singlets and triplets 
\{tma}) = ®{;}|si) «){„„} It^a) with an excitation energy 
given by the number of triplets N[{tma}]- At finite A the 
action of the ladder-Hamiltonian on these product states 
can be read off easily from (||): (i) H2 creates (destroys) 
pairs of nearest-neighbor (NN) triplets of equal a-index, 
(ii) given a pair of NN sites in a relative state of one 
singlet and one triplet H2 generates NN hopping of the 
triplet via exchange of the singlet and triplet, and finally 
(iii) H3 induces NN interactions between triplets. Con- 
tributions (ii) and (iii) define a three-flavored hard-core 
Bose system with NN hopping as well as NN interac- 
tions. These contributions renormalize the spectrum to 
0(A). Process (i) does not conserve the triplet number 
and changes the spectrum only to O(A^) at A <C 1. The 
essential physics of the bound states and optical absorp- 
tion is independent of 0(A^)-terms. They will be ne- 
glected hereafter. Based on this simplification it is con- 
venient to split the Hamiltonian into a bare 'kinetic' part 
Hi + H2 and a two-particle interaction H^. The bare 
one-triplet eigenstates \ka) of momentum k are 



\ka) 
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i/2|fca) - (|i?|) = l + Acos(fc) (4) 



where {k'a\kf3) — SapSkk' ■ This agrees with seminal work 
on spin-ladders P-H]. Hereafter the ground state energy 
{\H\) is set to zero. The bare two-triplet eigenstates in the 
singlet sector must incorporate the symmetry \tiatmf3) = 
{trnptia) ^s well as the hardcore constraint \tiatip) — 



\kq,S) 
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,ik(l+m)/2 



sgn(/ — m)sm(q(l ~ m))\lm, S)] 
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Ekq = 2[1 + Acos(fc/2) cos(g)] 



(6) 



where \lm,S) — '^^\tiat,na)/V3 refers to the singlet 
combination of two triplets in real space and Skq is 
the bare two-particle kinetic energy. k{2q) is the to- 
tal(relative) momentum, {k'q'\kq) — Skk'Sqq' and k G 
] — tTjTt], q g]0,7r]. (Anti)periodic boundary conditions 
apply to k{q) and \k,q) = —\k,—q). A discussion of the 
remaining two-particle states in the triplet and quintet 
sector is deferred to appendix W. 

After these preliminaries we are in a position to eval- 
uate the optical absorption due to phonon- assisted two- 
magnon emission (PME). This absorption results from 
a magnetoelastic excitation induced by the incoming 
photon-field and has been considered first in the planar 
geom etry o f the Cu02 square-lattice of the HTc mate- 
rials [|2 13|. Its microscopic justification can be applied 
with little change to the case of spin-ladders with corner- 
sharing Cu-0 structure. For the effective coupling one 
obtains 



Hls = -E{t) ^ [piu^oi +PA{'2Ufj,oi 



(7) 



E{t) is the time-dependent electric field which is polar- 
ized along the ladder. u^o(l,r)i denote displacement co- 
ordinates of the oxygen site on the /x-th leg of the lad- 
der and refer to the oxygen sites in cell 0{rigin) = I, 
L{eft) = I - 1, and R{ight) = 1 + 1. The effective 
charges pi and pA can be derived from an expansion 
up to first order in ui and E{t) of the Cu-O-Cu su- 
perexchange coupling. Expressions for pi and pA can 
be found in [n2lll3]. For cuprates which are strongly co- 
valent materials and for E{t) polarized parallel to the 
ladder, \pa\ ^ \pi\ is a reasonable approximation p3[ . 
Although the specific form of (0) depends on the corner- 
sharing Cu-0 structure it is conceivable that a similar 
PME is also possible in other spin-ladder compounds. 
Introducing the Fourier transformed phonon-operators 



Wfe/i 



^fc^ 



-k^ 



)/y/2Mu)k with the phonon disper- 
sion Wk and the reduced O mass M and replacing S'"j 
with the bond variables Sf and T" we get 



^is = i^W^ EKp + «U^)^fc 



V2 
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have 



where terms proportional to S^S"_^_-^ and S"Ti°^^ 
been dropped since they do not act on a pure singlet 
product-state. Photon momentum conservation has been 
neglected since the wave length of light is large compared 
to the lattice spacing. 
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FIG. 1. Dispersion of the S = 0,1, and 2 bound states 
relative to the bare two-triplet continuum. 

The zero-temperature absorption coefficient I{uj) at 
frequency u, normalized to the electric field intensity 
jE'P, is obtained from Fermi's golden rule. To simplify 
technical matters we assume only optical phonons of fre- 
quency ujk — ^0 to be involved in the absorption pro- 
cess. Since our model discards any intrinsic magnetoe- 
lastic coupling the phonons represent a mere 'momentum 
bath' allowing for arbitrary recoil of the two triplets |14| . 
Tracing over the phonons leads to 
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2TT}^\{fk\Pk\0)\^5iLj - Efk - u;o) 
-2Im J2mPl\kq,S) 
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where z = w + zO^ — cjo and for the remainder of this 
work we redefine oj — uji^^uj. While \fk) are the dressed 
two-particle eigenstates of the ladder of total momentum 
k, spin zero and energy Efk, \kq, S) are the bare states 

Next we have to evaluate the resolvent (^_1|). Since 
\kq,S) contains only two triplets and, as noted previ- 
ously, we neglect the coupling between sectors of different 
triplet number this task is equivalent to solving an inter- 
acting two-particle problem which can be done exactly. 
To this end we observe that 

{kq,S\H^\k'q',S) = -^^4fc'Asin(g)sin(g') (12) 

{k'q, 5|P, |0) = -./IpTpfce-^^'/s sin(<z)4fc' (13) 
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This can be combined with (O) to sum the T-Matrix 
series for {kq, S\l/{z — H)\kq' , S) yielding 
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FIG. 2. Optical absorption spectrum for pi = 
identical spectrum, y-axis stretched by 100. 
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where the bare two particle resolvent G{k, z) is given by 



G{k,z) 



dq 



sin^(q) 



ATrJ_^ "z- (2 + 2Acos(fc/2)cos(g)) 



sign[Re(a)]A/c 



1 



4Acos(fc/2) 



(15) 



with a = (2 — z)/[2Acos(fc/2)]. The zeros at a; = E^ and 
momentum k of the denominator on the r.h.s. of ( |l4| ) 

1 + 4AG(fc, £;f ) = (16) 

correspond to (anti)bound states in the two-triplet spin- 
singlet sector. Since both, the boundaries of the two- 
particle continuum (||) as well as the solutions of ( p^ ) 
are functions of k only if energies lu are rescaled in terms 
of the quantity (w — 2)/A, the fc-space structure of the 
rescaled two-particle spectrum is independent of A within 
our approximation. In particular, a singlet bound-state 
exists at all k in the Brillouin zone. At fc = tt its binding 
energy is largest with E^ = 2 — A while at fc = the 
binding energy is zero. This is shown in fig. nl which 
depicts the bound-state dispersion along with the two- 
particle continuum. Also shown are triplet (S — I) and 
quintet (S = 2) (anti)bound-states, which are optically 
inactive. These states are discussed in appendix ^. Re- 
cently, similar collective states have been predicted using 
bond-boson techniques |15|-[l7t]. 

From (13 l^ , p^ ) we can obtain the absorption by nu- 
merical integration. The rescaled intensity XI{uj)/pa is 
a function of lo — {uj — 2)/A and pi /pa only and does 
not dependent on A. Figure y depicts the absorption 
spectrum for pi — Q. We note that our choice of the 
frequency variable Cj implies that zero incoming photon 
energy corresponds to the point — 2/A on the x-axis of 
fig. 0. For any finite spin gap, i.e. A < 1, this point is 
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FIG. 3. Optical absorption spectrum for pi — 
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off the range plotted which refers only to the frequency 
window —2 < (D < 2 in which absorption occurs. Figure 
g demonstrates that the bound state has a profound im- 
pact on the absorption spectrum which comprises almost 
completely of a structure due to the integrated density 
of states of the bound state. In particular there is a van- 
Hove singularity at w = 2 — A which results from the 
maximum in the bound-state dispersion at k — n. Only 
in the special case of p/ — — 4p^ this van- Hove singularity 
disappears. The inset in fig. focuses on the remaining 
spectral weight. It has a maximum at the center of the 
two-triplet continuum, however its weight is very small 
as compared to the contribution from the bound states. 
At finite values of pj a step-like feature appears at the 
spin-gap and for p/ < a dip occurs due to interference 
in the effective charge pk- At present we believe that 
such features are of academic interest only. Yet, for com- 
pleteness a typical spectrum displaying the latter effects 
is shown in fig. 0. 

The preceding discussion does not rule out an inter- 
pretation of van-Hove like structures in an experimen- 
tally determined absorption spectrum solely in terms the 
momentum integrated bare two-triplet continuum rather 
than in terms of bound states. However, taking this point 
of view and using that \pi/pa\ ^ 1, which is most likely 
the case, the absorption spectrum would rather be sym- 
metric with respect to the center of gravity of the mea- 
sured spectrum. Finally we note, that in case of a finite 
dispersion of the phonon additional structures may ap- 
pear which are not contained in our present analysis. 

In conclusion we propose that singlet bound-states 
have a dramatic impact on the optical absorption of 
dimerized spin-ladder systems which allow for phonon- 
assisted multi-magnon absorption of the Lorenzana- 
Sawatzky type. We hope that our work may stimulate 
further experimental investigation of this issue. 
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APPENDIX A: S=l,2 STATES 

The triplet and quintet two-particle states in real space 
correspond to the combinations \lm,a) — J^Bj £q/37 

\tlf3tm'y)/V^ and \lm,aP) = \^^W^{\tlatm0) + \tiptma) " 

(5q./3(2/3)^ |ti^im7)) respectively. Due to spin rota- 
tional invariance we may focus on the total-S*^ = 
components only, i.e. \lm, P) — \lm,z) (triplet) and 
\lm,D) = \lm,zz) (quintet). The bare two-particle 
eigenstates are 

\kq, P) = ^J2 e*['=('+'"'/2+«('-™)l \lm, P) (Al) 



\kq,D) = , ^ :Y[i 



,ik{l+m)/2 



sgn{l — m)sm{q{l — m))\lm,D)] (A2) 

where identical properties regarding the momenta as in 
(h) apply and the bare two-particle energies of the states 
\kq,P) and \kq,D) are identical to Skq of ^. Note that 
both, P- and Z?-states conform with the hardcore con- 
straint \tiatip) = 0. From the identity (S;-HS„i)^/2 — 2 — 
S/ • Sm and from (n2|) it follows directly that 



(kq,X\H3\k'q',X) 



cx 



N - 1 



5kk'>^sm{q)sm{q') (A3) 



with Cx = 4, 2, and —2 for X = S, P, and D respectively. 
Therefore the energies E^ of the (anti) bound states in 
all spin channels are obtained from the zeros of the single 
equation 



l + cx\G{k,E^) = 



(A4) 



with G{k, z) as in (|l5|) and cx set according to the total 
spin. For X = P{D) a(n) (anti)bound state exists for 
k > 27r/3 for all values of A with E^^"^ = 2-{+)X/2. In 
fig. W the dispersion of these (anti)bound states relative 
to the two-triplet continuum is depicted. 
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